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probably  a  consequence  of  the  seasonal  intrusion  into  the  region  of  the  low  salinity  plume  of  the  Yangtze  River.  Monitoring  changes  in  the  dispersion 
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This  article  p  resents  find  tugs  regarding  the  origination  and  evolution  of  a  large  train  of  internal  solitary  waves  as 
observed  through  Synthetic  Aperture  Radar  (SAR)  The  internal  wave  train  understudy  propagates  into  deeper  waters 
away  from  the  coast  and,  most  important,  appears  to  he  generated  at  a  well-defined  point  along  a  steep  topographic 
variation.  The  shallow  bottom  topography  revealed  itself  in  theSAR  image  and  provided  the  opportunity  to  determine, 
with  a  large  level  of  accuracy  ,  the  precise  area  where  an  internal  bore  originated.  The  bore  traveled  east  and  broke 
up  into  a  train  of  solitary  waves  The  packet's  phase  speed,  wave  amplitude,  mixed-layer  depth,  and  density  difference 
for  a  two-layer  structure  were  remotely  estimated  based  on  measurements  made  from  the  SAF  image,  the  knowledge 
of  the  local  bathymetry,  historical  climatological  data,  and  a  two-layer  mode)  Farther  along  the  propagation  path  the 
internal  waves  exhibited  an  extremely  large  degree  of  dispersion,  which  wot  related  to  changing  stratification  con 
ditions  along  the  travel  direction  that  are  probably  a  consequence  of  the  seasonal  intrusion  into  the  region  of  the  low 
salinity  plume  of  the  Yangtze  River  Monitoring  changes  in  the  dispersion  characteristic#  of  a  aoliton  train  from  one 
tidal  period  to  the  next  could,  therefore,  render  remote  information  regarding  changing  stratification  conditions  in 
the  area  This  particularly  long  train  of  internal  waves  consisted  of  seven  packets  that  traveled  Tor  three  days  along 
a  well-defined  direction.  A  new  analytical  model,  termed  the  Dnoidal  model,  was  used  to  describe  the  development  of 
the  internal  waves  with  time,  The  extended  life  of  this  phenomenon  provided  an  exceptional  opportunity  to  monitor 
and  model  iu  evulullon  as  well  as  to  assess  the  effects  of  dispersion  as  the  solitary  internal  waves  propagated  east, 
□  way  from  the  China  coast. 


ADDITIONAL  KEY  WORDS:  in  ter  not  wove  modch,  Synthetic  Aperture  Rodar  SAR . 


INTRODUCTION 

A  500-km  x  500-km  ( ScanS AR  wide  mode)  image  (Figure 
l)  was  used  to  monitor  the  origination,  evolution,  and  dissi¬ 
pation  of  a  large  train  of  internal  wave  packets  propagating 
seaward  off  the  China  coast.  The  images  were  acquired  Au¬ 
gust  22,  1998  at  0942  by  the  Canadian  satellite  Radarsal  1. 
Medium  resolution  images,  at  100  m  per  pixel,  and  full  res¬ 
olution.  at  50  m  per  pixel,  were  analyzed.  Seven  packets  of 
internal  bore  and  internal  wave  type  features  were  followed. 
The  internal  bores  were  generated  beside  a  conspicuous  to¬ 
pographic  feature.  The  internal  bores  that  disintegrated  into 
solitons  traveled  far  into  the  midseclion  of  the  Yellow  Sea. 
Observing  that  the  spacing  between  the  wave  packets  had  a 
high  degree  of  regularity  along  a  well-defined  direction  and 
that  the  overall  shape  of  the  wave  front  was  preserved  along 
the  travel  paLh,  corroborated  the  fact  that  this  was  actually 
a  train  of  internal  waves  originated  at  the  same  location.  Fur- 
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thermore,  the  effects  of  dispersion  were  dearly  seen  ns  the 
internal  waves  traveled  away  from  their  generation  point 

Unlike  previous  literature,  which  mostly  describes  internal 
waves  propagating  inshore  (Liu  et  al.t  1998;  Porter  and 
Thompson,  1998),  the  internal  wave  train  understudy  prop¬ 
agated  eastward  into  deeper  waters  away  from  the  coast  and, 
most  important,  appeared  to  be  generated  at  a  well-defined 
area  along  the  shelf  break.  The  remotely  inferred  profile  of 
bathymetry,  which  revealed  itself  in  the  SA.R  image,  was  con 
sistent  with  the  known  local  topography.  This  provided  the 
opportunity  to  determine  within  a  good  degree  of  accuracy, 
the  precise  area  where  these  internal  waves  originated.  The 
sleep  shelf  profile  seen  at  the  onset  of  the  internal  wave  train 
was  used  to  initialize  a  hydrodynamic  model  that  reproduces 
the  internal  bore  seen  dose  and  offshore  of  the  shelf  break 
(Warn -Yarn ah  ct  ait  2005),  We  are  reporting  hero  on  the 
results  obtained  with  a  newly  introduced  analytical  model 
(Apel,  2003),  which  is  used  in  this  work  to  help  assess  the 
evolution  of  internal  waves  as  they  traveled  east. 

The  packet's  phase  speed,  wave  amplitude,  wavelength, 
and  water  column  properties  were  studied  in  relation  to  mea¬ 
surements  contained  in  the  Synthetic  Aperture  Radar  (SAR) 
image  and  with  the  knowledge  of  the  local  bathymetry  and 
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Figure  1.  A  section  of  the  Yellow  Sea.  depicting  the  focii*  aren,  ns  seen  in  a  IWarsotl  ScanSAH  hvirini  500  x  500  km  image.  Pixel  npudng  i*  100  m, 
August  22,  0040.  'fho  Inrgv  Sfimirimilftr  pattern  seen  M  the  western  pari  of  the  image  f  Figure  lb)  is  h  manifestation  of  the  underlying  bathymetry . 

Internal  wave  packets  arc  formed  at  this  steep  topography  und  travel  enst  The  Op  of  the  Korean  peninsula  and  Chgju  Inland  cun  he  seen  at  the  right, 
Figures  Id  ig  shown  the  signature  of  an  internal  boie-  as  it  leaves  the  shelf  break,  fUs  integrates  into  solitons,  mid  propagates  southeusL 


historical  dimaloingic&l  data.  The  signature  and  long  life  of 
this  event  provided  a  good  opportunity  to  monitor  and  model 
its  evolution  and  Lo  assess  the  effect  of  dispersion  as  the  in¬ 
ternal  waves  traveled  away  from  the  China  coast  (Figure  1). 


The  article  is  organised  by  first  providing  some  background 
regarding  SAR  studies  of  nonlinear  internal  waves  and  re¬ 
viewing  the  physical  settings  under  which  the  phenomenon 
takes  place.  The  overall  number  of  wave  packets  and  travel 


Jnumnl  of  Coufttfd  Research.  Vol.  22,  No  G.  20DB 


Internal  Solitary  Waves  in  the  Yellow  Sea 


1405 


direction  is  summarized.  The  topographic  feature  where  the 
first  interna)  bore  is  generated  is  studied  in  detail  and  its 
surface  signature  in  the  5AR  image  is  analyzed.  Phase  speed 
calculations  are  performed  assuming  forcing  by  the  semidi¬ 
urnal  tide.  Estimation  of  the  internal  wave  amplitude  is  ob¬ 
tained  from  the  soliton  half  width,  which  in  turn  is  derived 
from  the  distance  between  the  dark  and  bright  patterns  of  its 
surface  signature  in  the  SAR  image.  Based  on  the  local  ba¬ 
thymetry,  wave  phase  speed,  and  mixed-layer  thickness,  as 
determined  from  SAR  imagery,  an  estimation  of  the  density 
contrast  between  the  upper  and  lower  layers  is  obtained.  A 
two-dimensional  spectrum  of  the  SAR  image  permits  the 
quantification  of  the  packet  wavelength,  travel  direction,  and 
dispersion  effects.  The  large  degree  of  dispersion  seen  after 
the  packets  propagated  farther  south  is  associated  to  chang¬ 
ing  stratification  conditions,  possibly  related  to  the  seasonal 
intrusion  of  the  low-salinity  plume  of  the  Yangtze  River.  Fi¬ 
nally,  the  newly  introduced  Dnoidal  model  (Apel,  2003)  is 
used  to  study  the  evolution  of  the  internal  waves  as  they 
propagate  away  from  the  coast. 

INTERNAL  WAVE  DETECTION  USING  SAR 

It  is  known  that  tidal  flow  over  steep  topography  depresses 
the  stratification  interface  generating  an  internal  bore.  The 
bore  starts  propagating,  and  its  leading  edge  becomes  steeper 
through  nonlinear  effects.  With  time,  frequency,  and  ampli¬ 
tude,  dispersion  set  in  and  the  leading  edge  of  the  bore  dis¬ 
integrates  into  solitons. 

Synthetic  Aperture  Radar  detects  the  variations  in  sea  sur¬ 
face  roughness  associated  with  internal  waves.  These  chang¬ 
es  in  surface  roughness  due  to  convergence  and  divergence 
translate  into  bright  and  dark  patterns  in  the  pixel  intensity 
profile  of  the  SAR  image.  For  depression  solitary  waves,  the 
bright  bands  in  the  SAR  signature  are  in  the  propagation 
direction,  Solitary  waves  that  elevate  rather  than  depress  the 
interface  show  dark  bands  in  the  propagation  direction  (Liu 
et.  aLt  19981 

According  to  the  experiments  carried  out  by  Ramirez  and 
Renouard  (1998)  in  a  13-m  tank,  the  ebbing  tidal  currents 
are  responsible  for  the  formation  of  internal  waves  at  the 
shelf  break.  These  waves  are  observed  propagating  inshore, 
over  the  shelf,  and  offshore  away  from  the  shelf.  The  internal 
waves  under  study  offer  a  remarkable  opportunity  to  confirm 
that  the  lee-wave  formed  by  a  strong  ebbing  flow  over  topog¬ 
raphy  can  actually  lead  to  offshore  propagating  waves.  A 
model  by  Lamb  (1994)  illustrates  the  generating  mechanism. 

REGIONAL  CHARACTERISTICS 

The  Yellow  Sea  has  an  average  depth  of  44  m  and  a  max¬ 
imum  depth  of  about  100  m*  In  the  summer,  when  the  strat¬ 
ification  is  strongest,  the  mixed  layer  shallows  to  the  range 
of  10  m  to  20  m.  SAR  images  acquired  during  that  season 
show  a  myriad  of  internal  wave  signatures  approximately  be¬ 
tween  33°  N  to  36°  N.  The  focus  area  extends  from  33°12'  N 
to  34°12'  N  and  from  122*6'  W  to  125*30'  W,  The  soliton  train 
under  study  propagates  from  about  34°  N,  122°  W  to  33*24' 
N,  124*30'  W, 


The  Radarsat  data  set  parameters  were  the  following: 

•  Date:  August  22,  1998 

•  Time:  0942:25 

•  Orbit  Number:  14604 

•  Satellite  ID:  Radarsat- 1 

•  Processing  Site:  Alaska  SAR  Facility 

•  Size:  57,344  kB 

•  Beam  Mode:  ScanSar  Wide — 500  km  X  500  km 

•  Beam  Position:  SWB 

•  Product  Type:  Geocoded 

•  Ascending/Descending:  Ascending 

•  Pixel  Spacing:  100.00  m 

•  Northernmost  Latitude:  36°8'  16.44" 

•  Southernmost  Latitude:  31*3*13.32” 

•  Westernmost  Longitude:  1 21*13 '4.8* 

•  Easternmost  Longitude:  127*2 '3.48* 

During  August,  the  low  salinity  plume  from  the  Yangtze 
river  can  intrude  into  this  region  and  cause  an  increase  in 
vertical  stratification  (Hur,  Jacobs,  and  Teague,  1999).  The 
Yellow  Sea  Cold  Bottom  Water  is  another  water  mass  that 
can  also  intrude  into  this  region  as  it  moves  southward 
(Park,  1986).  Therefore,  the  trains  of  internal  solitary  waves, 
as  they  propagate  southwest,  will  not  only  encounter  changes 
in  depth  but  also  are  likely  to  come  across  changing  strati¬ 
fication  conditions  as  well. 

INTERNAL  BORE  GENERATION  AND 
ASSOCIATED  BOTTOM  TOPOGRAPHY  AS  SEEN 
THROUGH  SAR 

The  possibility  to  “see'’  shallow  topographic  features  at  the 
bottom  using  SAR  images  is  achievable  because  strong  tidal 
currents  flowing  over  bottom  topography  cause  a  modulation 
in  the  surface  current  velocity.  This,  in  turn,  causes  varia¬ 
tions  in  the  spectrum  of  wind-generated  waves,  which  show 
up  subsequently  as  intensity  modulations  in  the  SAR  images 
(Hsu,  Liu,  and  Mitnik,  2001), 

The  M2+S2  tidal  amplitudes  at  the  shelf  break,  where  the 
internal  bore  under  study  is  formed,  range  between  1.5  m  to 
2  m  whereas  the  tidal  currents  range  between  0.6  m/s  and 
1.2  m/s.  The  flooding  and  ebbing  tide  flowing  over  a  steep 
topographic  feature  depresses  the  stratification  creating  an 
internal  bore.  Note  that  the  shape  of  the  bore  follows,  to  a 
certain  extent,  the  shape  of  the  shelf  break  at  which  it  was 
generated  (Figure  lc)*  The  underwater  topography  reveals  it¬ 
self  through  quasistationaiy  patterns  of  bright  and  dark  pix¬ 
els  bands  (looking  from  west  to  east).  The  bore,  instead, 
shows  an  opposite  pattern  of  intensity  that  goes  from  dark  to 
bright.  Because  the  bore  is  traveling  east,  it  should  be  de* 
pressing  the  stratification.  As  mentioned  in  Lru  et  at .  (1998), 
a  disturbance  of  depression  can  yield  a  pattern  of  bright  pix¬ 
els  in  the  direction  of  travel.  In  turn,  the  shallow  underwater 
bathymetry  and  the  associated  tidal  flow  seem  to  be  intro¬ 
ducing  a  disturbance  of  elevation  in  the  stratification  that 
gives  the  exact  opposite  pattern  of  convergence-divergence 
at  the  sea  surface.  This  translates  into  opposite  patterns  of 
pixel  intensity  in  the  SAR  image.  Further  information  re¬ 
garding  the  theory  of  radar  imaging  of  internal  waves  and 
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Figure  2.  Note  (hi*  relatively  *U*'p  Kiope  anumd  24  'N.  122' IS'  E,  jiwl  where  the  intermit  bore  is  generated  The  bnlhymdi)’  profile  n(  the  mIc  could 
related  le  what  appears  to  be  a  refraction  of  the  in  ten  ml  bore  shortly  after  its  formation,  Depth  scale  is  in  meters. 


surface  signature  of  bottom  topography,  including  their  re¬ 
lation  to  wind  and  look  direction,  can  be  found  in  Audits 
1 19851  and  Romkiskk  and  Am»kkk  ( 1997l.  The  appearance  of 
the  bathymetry  as  it  in  revealed  through  BAR  (Figures  lb  and 
Id  )  can  be  compared  with  the  actual  one-minute  bathymetry 
(Choi,  1999)  in  Figures  la  and  2,  A  three-dimensional  ren¬ 
dering  of  this  topographic  feature  was  also  generated  to  bet¬ 
ter  understand  the  phenomenon.  The  finger- like  structure 
seen  in  the  SAR  images  is  clearly  depicted  in  the  topographic 
variations,  A  typical  depth  profile  along  a  transect  perpen¬ 
dicular  Ln  the  shclT  break  is  illustrated  in  Figure  3. 

PACKET  CHARACTERISTICS  AND  PARAMETERS 

The  ScanSAR  image  shows  seven  packets  of  solitary  waves 
propagating  271  km  during  3  days  along  a  well-defined  di* 
rod  ion.  The  first  soli  ton  train  extends  along  its  crest  for 
about  132  km.  The  average  distance  between  packets  is  43 


km,  and  the  actual  location  of  each  packet  is  summarized  in 
Table  1.  The  average  phase  speed  is  1,01  m/s.  The  shape  of 
the  packet  wave  fronts  are  preserved  remarkably  well  from 
the  second  to  the  fifth  tidal  period. 

An  interesting  feature  of  the  particular  train  under  study 
is  that  after  one  semidiurnal  tidal  period  the  internal  bore, 
which  was  formed  near  the  shelf  break,  failed  to  break  up 
into  solitons,  Because  the  original  image  had  a  resolution  of 
100  m/pixel,  a  full -resolution  ScanSAR  image  with  a  resolu¬ 
tion  of  50  m/pixel  was  analyzed.  The  mentioned  internal  bore 
still  did  not  show  any  further  details  at  this  higher  resolu¬ 
tion,  The  fact  that  the  boro  appears  to  fail  in  breaking  up  into 
solitons  after  one  tidal  period  could  be  because  of  the  lack  of 
image  resolution  or  it  could  l>e  the  result  of  the  vary  ing  un¬ 
derlying  bottom  topography  near  the  generation  point  As  a 
matter  of  fact,  right  after  the  internal  bore  is  formed,  u  prop¬ 
agates  toward  a  sudden,  sill-type  topographic  feature  <Fig- 
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Figure  3.  Depth  profile  follows  the  central  propagation  path.  After  the 
shelf  break,  there  la  a  conspicuous  underwater  feature. 


uxes  2  and  3),  whose  interaction  with  the  bore  could  be  the 
subject  of  future  studies. 

Another  observation  is  that  the  second  bore  is  not  parallel 
to  the  first.  This  refraction  effect  could  be  related  bo  the 
northern  end  of  the  internal  bore  propagating  over  deeper 
waters  than  the  southern  end.  The  phase  speed  is  between 
0.8  m/s  and  1.03  m/s  depending  on  what  end  (south  or  north) 
is  measured  (Figure  4).  The  third  event  is  also  refracted  (Ta¬ 
ble  1)  with  respect  to  the  second.  After  this  initial  refraction, 
the  packets  preserve  their  overall  direction  for  the  following 
four  tidal  periods.  From  the  first  to  the  second  wave  packet, 
there  is  a  dear  increase  in  phase  speed.  Changes  in  the  shape 
of  the  packets  could  be  also  because  of  variations  of  the  ver¬ 
tical  stratification  along  the  propagation  path,  changes  of 
depth,  advection  by  a  variable  current  field,  and  variations 
of  the  tidal  forcing  from  one  cycle  to  the  next  (Porter  and 
Thompson,  1999). 

Along  the  travel  path,  there  is  a  depth  gradient,  from  shal¬ 
low  to  deep  water,  not  only  in  the  propagation  direction  but 
also  in  the  along-crest  direction  with  the  southern  end  of  the 
trains  being  over  shallower  waters  than  their  northern  ends. 
This  probably  accounts  for  the  fact  that  most  wavelengths, 


isw 


Figure  4,  Main  characteristicB  of  the  internal  wave  train. 


inferred  from  the  SAR  image,  are  smaller  toward  the  south 
end  of  the  train  fronts  than  toward  the  north  end.  In  addition, 
note  that  as  the  south  end  of  the  wave  crest  enters  shallow 
waters,  its  signature  disappears.  In  this  article,  we  will  show 
that  the  joint  consideration  of  packet's  structure,  propagation 
direction,  pattern  of  pixel  brightness  in  the  SAR  imagery,  and 
extent  of  dispersion  in  conjunction  with  any  observed  refrac¬ 
tion  permits  us  to  speculate  regarding  the  overall  profile  of 
bottom  topography,  variable  current  fields  that  might  be  ad- 
vecting  the  soliton  trains,  and  spatial  density  distribution  in 
the  region. 

MEASURING  INTERNAL  WAVE  PARAMETERS 
FROM  SAR  IMAGERY 

There  are  four  parameters  that  can  be  directly  measured 
from  SAR  images:  the  packet's  phase  speed,  the  soliton  half¬ 
width,  the  wavelength  at  the  trailing  ends  of  the  packets,  and 
the  stretched  wavelength.  The  phase  speed  can  be  calculated 
assuming  forcing  by  the  M2  tide  (12.42  h)  and  measuring  the 
distance  between  the  packets  fronts.  In  this  particular  case 
study,  there  is  little  variation  in  the  phase  speed  from  one 
packet  to  the  next  for  the  duration  of  the  different  tidal  pe¬ 
riods.  This  is  an  indication  that  the  background  current  speed 
is  fairly  constant  along  the  propagation  path. 

The  halfiwidth  length  scale  L  of  the  solitons  can  be  ob¬ 
tained  from  the  distance  that  exists  between  the  patterns  of 


Table  l,  Evolution  of  a  packet  of  internal  waves  as  it  propagates  during  three  days  in  the  Yeltcuv  Sea. 


Soliton  Train  # 

Bore  t 

Bore  2  Refracted 

Packet  l  Refracted 

Packet  2 

Packet  3 

Packet  4 

Packet  6 

Time  (h) 

0 

12,4 

24.B 

37.2 

49  6 

62 

74.4 

Distance  (km} 

36-46 

72-93 

132 

178 

226 

272 

Phase  speed:  C  (m/a> 

0.8-1.03 

0.S-1.05 

1.09B 

1,03 

1,075 

1,03 

Dominant  Wavelengths: 

A1  (m) 

530 

530 

530 

530 

A2  <m) 

630 

670 

670 

670 

A3  (m> 

930 

1000 

930 

A4  (m) 

1200 

1200 

A5  (ml 

1600 

1600 

A6  (ml 

2300 

2400 
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Figure  5.  Relation  between  the  iaopycnal  displace  me  lit  caused  by  the 
internal  solitary  wave  and  the  pixel  intensity  pattern  in  the  SAH  image. 
'Hie  distance  between  the  dark  and  bright  spots  is  termed  D  The  soli  ton 
half- width  length  scale  is  termed  L 


bright  and  dark  pixels  in  the  SAR  imagery.  These  pixel  ar¬ 
rangements  correspond  to  areas  of  peak  convergence  and  di¬ 
vergence  at  the  sea  surface.  Convergence  and  divergence  at 
the  sea  surface  are  related  to  the  surface  strain,  which  is  the 
horizontal  gradient  of  the  surface  velocity.  The  analytical 
functional  profile  of  the  surface  current,  in  turn,  has  the  same 
sech3  outline  as  the  isopycnal  displacement  (eee  Equations  1 
and  2)  caused  by  the  internal  wave  but  with  opposite  phase. 
The  peaks  and  the  troughs  in  the  surface  strain  give  the  peak 
in  convergence  and  divergence  and  can  be  found  by  using  the 
nulls  of  the  first  derivative  of  the  strain  rate,  which  is,  there¬ 
fore,  the  second  derivative  of  the  horizontal  current.  Hence, 
these  points  of  maximum  convergence  and  divergence  corre¬ 
spond  to  the  inflexion  points  of  the  surface  current  functional 
profile.  Because  this  current  has  the  same  "sech3  *  shape  as 
the  isopycnal  displacement,  they  also  correspond  to  the  in¬ 
flexion  points  in  the  soliton  functional  representation.  Be¬ 
cause  the  surface  strain  rate  can  be  considered,  on  a  first- 
order  approximation,  roughly  proportional  to  variations  in 
the  pixel  intensities  in  the  image  (Liu  et  a/.,  1998),  we  can 
define  D  as  the  distance  between  the  inflexion  points  in  the 
soliton  profile.  This  is  equivalent  to  saying  that  D  is  the  dis¬ 
tance  between  the  bright  and  dark  bands  in  the  pixel  inten¬ 
sity  patterns  of  the  SAR  image.  Because  the  half- width 
length  scale  L  of  the  soliton  is  a  parameter  in  the  argument 
of  its  sech3  functional  profile,  it  can  be  easily  related  to  the 
distance  between  the  inflexion  points  of  the  curve  through 
differentiation  of  the  analytical  solution  for  the  horizontal  ve¬ 
locity,  The  positions  of  peak  divergence  and  convergence  are 
obtained  by  setting  the  second  derivative  of  velocity  zero. 
This  relation  gives  L  =  D/1.32  (Small  et  a h  19991  Given  our 
system  look  direction.  Figure  5  depicta  the  relation  between 
the  isopycnal  displacement  caused  by  the  internal  solitary 
wave  and  the  pixel  intensity  distribution  in  the  SAR  image. 
Further  discussions  regarding  SAR  imaging  of  internal  waves 


in  association  with  wind  and  look  direction  can  be  found  in 
Romeiser  and  Alpers  (1997).  The  inflexion  points  in  the 
isopycnal  displacement  are  correlated  to  the  maximum  and 
minimum  pixel  brightness.  Note  that  D  is  larger  than  L  and 
that  ZL  is  a  good  indication  of  soliton  width.  In  our  case 
study,  we  estimated  D  at  150  m  to  200  m,  which  gives  an  L 
between  113  m  and  151  m.  We  will  show  later  that  the  soliton 
half-width  can  be  used  in  combination  with  upper  and  lower 
layer  depth  to  obtain  an  estimate  of  the  isopycnal  displace¬ 
ment 

The  wavelength  at  the  trailing  end  of  the  packet,  X<,  re¬ 
mains  constant  as  the  packets  evolve  with  time  and  it  is  re¬ 
lated  to  the  isopycnal  displacement  and  to  the  parameters  of 
dispersion  and  nonlinearity,  The  trailing  end  wavelength  was 
estimated  at  between  300  m  and  350  m.  The  distance  be¬ 
tween  successive  soli  tons  along  the  packet  is  termed 
^stretched  wavelength/*  X,  and  it  is  dependent  on  the  degree 
of  nonlinearity.  The  stretched  wavelength  is  always  larger 
than  the  wavelength  at  the  trailing  end  of  the  packet. 

MONITORING  THE  EFFECTS  OF  DISPERSION 

The  leading  wave  in  the  packet  has  the  largest  amplitude. 
Because  the  soliton  phase  speed  is  related  to  its  amplitude, 
the  leading  wave  in  the  packet  travels  faster  than  the  rest 
(Porter  and  Thompson,  1999).  Therefore,  the  distance  be¬ 
tween  a  given  wave  crest  and  the  next  decreases  from  the 
front  to  the  back  of  the  packet  Hence,  the  solitons  propagate 
along  the  density  interface  as  an  internal  wave  packet  with 
rank-ordered  amplitudes  and  wavelengths. 

A  two-dimensional  Fourier  spectrum  of  a  20-km  x  20-km 
section  of  the  SAR  image  was  taken  for  each  wave  train  along 
central  transects.  As  the  second  internal  bore  propagates 
away  from  the  shelf,  it  breaks  up  into  solitons  where  we  can 
see  a  tightly  packed  soliton  train  with  a  leading  wavelength 
in  the  order  of  530  m  to  670  m  (Figure  le).  The  surface  spec* 
tra  routine  gives  wavelength  in  meters  and  propagation  di¬ 
rection  with  180°  of  ambiguity  (Figure  6a).  The  surface  spec¬ 
tra  corresponding  to  the  packet  shown  in  Figure  If  is  depicted 
in  Figure  6b.  After  several  tidal  periods,  the  spectra  shows 
wavelengths  ranging  from  530  m  or  less  at  the  trailing  end 
of  the  packet  to  2400  m  at  the  front  end  (Figures  lg  and  6c). 

As  the  packet  propagates,  the  surface  spectra  dearly  shows 
the  effects  of  dispersion  with  distance  (time).  Table  1  sum¬ 
marizes  the  dominant  wavelengths  in  the  packets,  as  seen  in 
the  two-dimensional  Fourier  spectra  us.  distance  (time).  In 
general,  the  leading  wavelength  measured  directly  from  the 
SAR  image  is  larger  than  the  largest  detected  by  the  two- 
dimensional  spectra.  This  is  because  the  image  is  coarse,  and 
sometimes  several  secondary  wavelengths  measure  about  the 
same,  thus,  concentrating  the  most  spectral  energy  at  the 
particular  wavelength.  In  this  regard  Levine  and  Boyp 
( 1998)  noted  that,  in  practice,  we  can  usually  find  a  few  boI- 
itons  within  a  packet  that  are  not  rank  ordered. 

For  the  first  three  events,  the  effects  of  refraction  were 
pointed  out.  For  the  last  four  packets,  measures  were  taken 
along  a  representative  transect  close  to  the  center  of  the  prop¬ 
agation  path.  The  wavelengths  at  the  trailing  end  of  the  pack* 
ets  tend  to  remain  the  same  whereas  the  wavelengths  near 
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Figure  <»,  The  wavelengths  m  t  he  first  train  urr*  tightly  packet  nml  dhow  little  sign  of  dispersion,  This  cart  Ik*  seen  in  the  actual  rudnr  imagery  <  Figure 
ho  ns  well  ns  in  the  surface  wave  spocim  of  the  soli  ton  trtnn  tuh  Wavelength  is  given,  in  meters,  by  the  concentric  circles  Figures  Jbl  and  i  c  i  depict  the 
surface  wave  spccim  showing  the  effects  of  diapersinn  m  the  pockets  shown  at  Figure  If  and  Figure  lg. 


the  leading  edge  tend  to  increase'.  This  is  because  the  centroid 
of  the  packet  travels  at  the  long-wave,  linear  phase  speed  tv, 
whereas  the  front  of  the  packet  travels  at  the  nonlinear  phase 
speed  c  <APKh,  2003 1.  The  nonlinear  phase  speed  is  larger 
than  the  linear  phase  speed,  and  this  accounts  for  the  in- 


\ 

creasingly  larger  wavelengths  observed  at  the  leading  edge 
of  the  packets. 

The  amount  of  dispersion  observed  in  these  trains'is  sig¬ 
ns  fkunL  and  because  different  trains  at  difTereni  times  show* 
different  degrees  of  dispersion,  it  is  possible  to  conclude  that 
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changing  conditions  in  stratification  and  depth  or  advection 
by  variable  current  fields  combined  with  the  nonlinear  char¬ 
acteristics  of  the  phenomenon  could  be  assisting  the  extreme¬ 
ly  large  spread  in  wavelength  seen  at  the  front.  Note  that 
from  the  second  to  the  third  packed  the  wavelengths  at  the 
leading  end  of  the  packet  almost  double  (Figures  If  and  Ig). 
This  is  consistent  with  a  suggestion  that  propagation  prop¬ 
erties  could  he  sensitive  to  slight  variations  in  the  stratifi¬ 
cation  (Levine  and  Boyd,  1998). 

ASSESSING  WATER  COLUMN  PROPERTIES  RASED 
ON  SAR  IMAGES  OF  INTERNAL  WAVES  AND  A 
TWO-LAYER  MODEL 

Given  the  internal  wave  characteristics  described  above, 
the  known  bottom  topography,  and  the  historical  climatolog¬ 
ical  data,  changes  in  vertical  stratification  along  the  propa¬ 
gation  path  could  be  remotely  inferred  from  the  SAR  data. 
These  could  include  variations  in  the  mixed-layer  depth  as 
well  as  changes  in  the  fractional  density  difference  between 
the  layers.  Although  the  stratification  in  the  summer  is 
strong  with  the  mixed  layer  at  a  relatively  constant  depth, 
as  the  train  of  solitary  waves  travel  for  3  days  along  a  272- 
km  path,  they  are  likely  to  encounter  variable  depth  and 
stratification  conditions.  This  particular  propagation  path  is 
dose  to  the  boundary  of  the  Yellow  Sea  and  the  northern  East 
China  Sea.  During  the  month  of  August,  stratification  con¬ 
ditions  could  vary  along  this  transect  due  to  uneven  intrusion 
of  different  water  masses  such  as  the  low-salinity  plume  from 
the  Yangtze  River  combined  with  the  Taiwan  Warm  Current 
Water  (Yangtze  River  Diluted  Water)  and  the  Yellow  Sea 
Cold  Bottom  Water  (Park,  1986;  Su  and  Wenc,  1994). 

Inference  of  Mixed  Layer  Depth  Behavior  From  SAR 
Imagery 

One  of  the  methods  reported  in  the  literature  (Porter  and 
Thompson,  1999;  Porter  et  ai ,  2001)  for  determining  the 
mixed-layer  depth  is  based  on  the  location,  as  inferred  from 
SAR  observations  of  the  region,  where  the  internal  wave  sig¬ 
nature  disappears.  This  technique  is  based  on  different  stud¬ 
ies  (Helfrich,  1992;  Helfrich  and  Melville,  1986;  Saf- 
farinia  and  Kao,  1996),  which  show  that  dissipation  of  in¬ 
ternal  waves  propagating  into  shoaling  waters  occurs  when 
both  layers  achieve  equal  widths.  Thus,  knowledge  of  the  ba- 
thymetry  at  the  location  where  the  internal  wave  signature 
can  no  longer  be  seen  allows  an  estimation  of  the  mixed* layer 
depth. 

Another  technique  can  be  inferred  from  the  processes  stud¬ 
ies  reported  by  Liu  et  at.  ( 1998),  in  which  depression  solitons 
that  are  propagating  into  shoaling  waters  change  into  soli- 
tons  of  elevation.  This  phenomenon,  where  the  solitons 
change  from  depression  to  elevation  disturbances,  occurs  af¬ 
ter  passing  a  location  at  which  both  layers,  of  a  two-layer 
ocean,  are  of  equal  width.  After  this  inversion  has  occurred, 
the  patterns  of  surface  divergence  and  convergence  at  the  sea 
surface  also  change,  and  the  pixels  in  the  propagation  direc¬ 
tion  become  dark  rather  than  bright.  SAR  observations  and 
numerical  simulations  of  this  process  can  be  found  in  Liu  et 
at. .  (1998), 


It  is  evident,  looking  at  the  BAR  image  (Figures  lb  and  Id) 
in  correlation  with  the  underlying  bathymetry  (Figure  2), 
that  the  internal  bore  has  originated  along  the  shelf  break 
starting  at  depths  that  are  in  the  range  of  30  m  to  40  m.  The 
shelf  edge  surface  signature,  however,  continues  south  over 
shallower  waters.  Because  there  is  no  apparent  reason  for 
this  to  happen,  other  than  the  existing  depth  and  stratifica¬ 
tion  conditions,  we  hypothesize  that  the  mixed-layer  depth 
could  be  near  the  critical  point  at  this  location.  This  allows  a 
rough  estimation,  which  should  be  half  the  depth  at  which 
the  disturbances  are  no  longer  seen,  of  approximately  20  m. 
This  speculation,  which  loosely  follows  the  first  technique  for 
mixed-layer  depth  determination,  is  in  very  good  agreement 
with  historical  climatology  (Navoceano,  1997)  for  the  re¬ 
gion. 

It  could  be  argued  that  the  surface  signature  of  an  internal 
wave  could  also  become  negligible  or  nonexistent  due  to  ad¬ 
verse  wind  conditions  and  not  necessarily  because  of  its  dis¬ 
integration.  It  is  highly  unlikely,  however,  that  this  sharp 
cutoff  in  the  surface  signature  seen  in  the  a  long-crest  direc¬ 
tion,  throughout  all  the  packets,  could  be  because  of  the  wind. 
If  this  were  the  case,  the  surface  signature  of  the  underwater 
topography,  which  is  roughly  parallel  to  the  trains,  should 
also  show  a  sharp  cutoff  in  that  direction. 

To  apply  the  second  methodology,  we  need  to  detect  a  de¬ 
finitive  change  in  the  pattern  of  pixel  brightness  along  the 
packet's  propagation  direction.  Although  no  strong  indication 
of  such  a  phenomenon  has  been  detected  in  the  image  under 
study,  there  is  an  internal  bore  in  the  proximity  of  the  40-m 
isobath  that  warrants  further  studies  through  the  use  of 
more  sophisticated  image- processing  techniques. 

A  Simplified  Two-Layer  Model 

Shallow  sea  internal  solitons  are  usually  modeled  based  on 
the  classical  Korteweg-de Vries  (K-dV)  equation: 

2 + (c# + aT,)S +  = 0  u) 

where  is  the  displacement  of  the  isopycnal  surface,  c0  is  the 
linear  phase  speed,  a  is  the  nonlinear  coefficient,  and  y  is  the 
dispersion  coefficient  The  K-dV  equation  can  be  solved  for  a 
two-layer  mode!  representation  of  the  ocean.  The  solution 
consists  of  an  invariant  squared  hyperbolic  secant  pulse  and 
a  Cnoidal-wave  represented  by  a  periodic  Cnoidal  Jacobian 
elliptic  function  (Apel,  2003;  Gerkema,  1994),  The  single 
soliton  form -preserving  solution  in  the  ocean  has  been  rep¬ 
resented  by  the  hyperbolic  secant  pulse  (Osborne  and 
Burch,  1980}  and  is  given  by  Equation  (2). 

(2) 

where  nc  ia  the  amplitude  or  peak  pynocline  displacement,  c 
the  phase  speed,  and  L  is  the  soliton  characteristic  half-width 
length  scale.  Equations  (3)  to  (8)  give  the  coefficients  of  the 
K-dV  equation  for  a  two-layer  model  with  upper  and  lower 
layer  widths  h  1  and  k2,  respectively  (Osborne  and  Burch, 
1980), 
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Table  2,  Summary  of  change*  m  depth,  phase  speed,  amplitude t  and  density  contrast  as  the  internal  iooum  propagate  coil  southeast  along  a  central 
transect. 


D  “  150  m 

L  *  H3.fi  m 

C  tm/aj  from  BAR 

Al  (m)  from  SAR 

Depth  (m)  from 
Bathymetry 
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20 

90 
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0.925 

20 

50 

3,7 

0.0068 

-0,0224 

89,7 

Packet  2 

1.098 

20 

60 

33 

0.0085 

-0.039 
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Packet  3 

1.03 

20 

60 

3.3 

0.0075 

-0.037 

131.89 

Packet  4 

1,075 

20 

60 

3.3 

0,0081 

-  0,0378 

137.6 

Packet  5 

1.03 

20 
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-0.044 
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Determination  of  the  Density  Contrast  between  Layers 
Based  on  SAR  Imagery  and  the  Two-Layer  Model 
Equations 

Using  the  half-width  of  the  solitons,  as  measured  from  the 
SAB  image,  assuming  that  the  upper  layer  roughly  preserves 
its  width  along  the  propagation  path,  and  considering  the 
bathymetry  for  the  region,  we  can  estimate  the  density  con¬ 
trast  between  both  layers.  Although  only  historical  data  are 
available  to  validate  our  results,  to  establish  the  usefulness 
of  this  technique,  the  two-layer  model  equations  and  obser¬ 
vations  made  from  SAR  will  be  used  to  infer  fractional  den¬ 
sity  difference,  (p*  -  PiVPa-  between  both  layers  along  the 
propagation  path. 

The  upper  layer  thickness  is  Al,  and  the  lower  layer  thick¬ 
ness  is  A 2  with  associated  densities  At  and  Aj,  the  density 
difference  is  Ap  =  pa  -  plr  and  the  fractional  density  differ¬ 
ence  or  normalized  density  is  Ap norm  =  (pa  -  p^l/p*.  From 
the  knowledge  of  the  phase  speed  c,  the  distance  between 
bright  and  dark  pixels  bands  in  the  SAR  image  D ,  and  the 
layer  thickness,  the  equations  describing  a  simplified  two- 
layer  model  can  be  used  to  calculate  the  density  contrast  be¬ 
tween  a  two-layer  structure. 

The  phase  speed  of  the  solitons  can  be  obtained  by  assum¬ 
ing  that  the  spacing  between  the  solilon  packets  is  the  result 
of  the  semidiurnal  tidal  forcing  (12.4  h)  and  by  measuring 
the  distance  between  the  wave  fronts.  We  have  done  this  in 
a  central  transect  along  the  propagation  path.  In  the  along- 
crest  direction,  there  is  a  steep  gradient  in  bathymetry. 

Calculations  are  illustrated  for  packet  number  three.  The 
third  packet  has  an  estimated  phase  speed  of  c  =  1.03  m/s 
along  a  central  transect,  the  distance  between  dark  and 


bright  pixels  was  estimated  at  D  —  150  m  from  SAR  imagery, 
the  bottom,  Al  +  A2,  was  estimated  at  60  m  from  the  ba¬ 
thymetry,  and  the  mixed -layer  depth  was  estimated  at  Al  «= 
20  m,  as  explained  in  a  previous  section.  Equations  (3)  to  (7) 
applied  to  the  third  internal  wave  packet  give  the  following 
results: 

L  =  113.6  m,  Tjo  =  3,3,  Apnorm  -  0.0075, 

c0  =  0*989,  a  «  -0,037,  y  -  131.89,  (8) 

The  overall  results  are  summarized  in  Table  2  for  the  dif¬ 
ferent  depths  and  phase  speeds  along  the  propagation  path 
with  the  mixed-layer  depth  assumed  to  be  constant.  As  the 
internal  waves  propagate  east  southeast,  the  density  contrast 
tends  to  increase  (Table  2),  and  because  there  is  no  large 
variations  in  depth  along  this  transect,  the  increase  could  be 
because  of  the  uneven  intrusion  of  different  water  masses 
along  the  propagation  path.  This  is  consistent  with  the  strat¬ 
ification  conditions  that  are  known  to  exist  during  the  month 
of  August  when  the  low-salinity  plume  of  the  Yangtze  River 
is  farthest  north  (Him,  Jacobs,  and  Teague,  1999). 

Note  that  from  the  second  to  the  third  packet,  the  leading 
wavelengths  almost  double  (Figures  If  and  lg).  This  could  be 
a  manifestation  of  the  highly  nonlinear  characteristics  of  this 
phenomenon  coupled  to  changing  stratification  conditions. 
An  increase  in  the  density  contrast  should  cause  an  increase 
in  the  nonlinear  phase  speed  at  the  leading  end  of  the  packet, 
and  as  a  consequence,  the  increase  in  wavelength  at  the  front 
could  be  larger  than  that  expected  under  unchanged  strati¬ 
fication  conditions.  Other  packets  arriving  to  the  area  from 
different  points  also  show  a  similar  large  degree  of  disper¬ 
sion,  thus  supporting  the  hypothesis  that  there  could  be  an 
intrusion  of  lighter  water  into  the  region. 

From  carefully  measured  phase  speed  and  depth  distribu¬ 
tions  at  various  points  along  the  internal  wave  train,  a  two- 
dimensional  density  contrast  field  of  the  area  could  be  con¬ 
structed.  When  the  lower  layer  density  is  known  from  histor¬ 
ical  climatology,  then  the  upper  layer  density  can  also  be  es¬ 
timated.  The  lower  layer  is  usually  more  stable,  and  its 
density  does  not  change  as  much  as  it  does  in  the  upper  layer. 
Even  though  the  absolute  values  of  the  estimates  for  the  den¬ 
sity  contrast  in  Table  2  appear  high,  it  is  likely  that  from  a 
comparative  point  of  view,  from  one  region  to  another,  or 
from  within  the  same  area  but  at  different  times,  they  could 
provide  useful  information  regarding  the  regional,  upper  lay¬ 
er,  relative  density  distribution. 

The  internal  wave  train  under  study  propagates  for  270  km 
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and  measures  up  to  130  km  in  the  along-crest  direction,  in* 
dicating  an  area  of  about  35,100  km*,  at  the  entrance  of  the 
Yellow  Sea,  which  could  be  remotely  monitored.  Note,  how¬ 
ever,  that  the  existence  of  this  particular  internal  wave  train 
could  be  seasonal  and  only  peak  spring  tides  might  support 
its  formation.  The  study  of  further  SAR  imagery  from  this 
area  is  needed  to  clarify  this  point  and  to  monitor  whether 
the  dispersion  extent  observed  throughout  the  internal  wave 
packets  is  related  to  seasonal  changes  in  stratification  con¬ 
ditions,  variable  current  fields,  or  other  regional  phenomena. 
Although  some  of  these  conclusions  are  speculative,  they 
serve  to  illustrate  the  potential  of  this  methodology  .  Coastal 
ocean  monitoring,  including  the  capability  to  identify  and  in¬ 
terpret  surface  features,  can  be  further  enhanced  when  SAR 
imagery  is  used  in  combination  with  optical  and  infrared  sen¬ 
sors  (JOHANNESSEN,  2000)* 

THE  DNOIDAL  MODEL 

The  recently  introduced  analytical  Dnoidal  model  (Apel, 
2003)  predicts  the  most  important  characteristics  of  conti¬ 
nental  shelf  solitons  as  they  evolve  in  time  and  space*  The 
model  is  based  on  a  solution  of  the  K*dV  equation  for  the 
amplitude  and  the  Taylor-Goldstein  equation  for  the  vertical 
structure*  Among  the  properties  that  can  be  monitored  with 
this  model,  we  have  the  number  of  solitons,  linear  wave 
speed,  waveform,  and  attenuation.  Synthetic  Aperture  Radar 
images,  which  provide  the  means  to  measure  important  in* 
temal  wave  parameters,  such  as  wavelength  and  half-width, 
can  be  used  in  conjunction  with  historical  climatological  data 
and  numerical  or  analytical  models  to  further  our  under¬ 
standing  of  the  dynamics  of  solitons  in  the  ocean. 

Model  Background 

The  solution  to  Equation  (1)  given  by  the  hyperbolic  secant 
profile  in  Equation  (2)  does  not  allow  the  addition  of  new 
solitons  pulses  as  the  packet  propagates,  does  not  provide  for 
the  dispersion  effects  as  soliton  trains  evolve  with  time,  and 
does  not  model  the  depression  of  the  isopycnal  displacement 
by  the  internal  tide.  ApeL  (2003)  used  a  little-known  solution 
to  Equation  fl)  to  model  the  mentioned  effects.  Gurevich 
and  Pitaevskii  U973)  originally  used  this  solution  to  de¬ 
scribe  collisionless  shock  waves  in  plasma*  This  solution, 
termed  the  Dnoidal  solution,  is  multiplied  by  an  ad  hoc  term 
that  models  the  recovery  of  the  isopycnal  surface  to  its  equi¬ 
librium  state.  It  U  also  an  asymptotic  solution  to  the  K-dV 
equation  and,  therefore,  not  valid  near  the  internal  wave 
origination  point  (Apel,  2003).  As  time  passes  the  model  adds 
one  oscillation  per  buoyancy  period.  According  to  the  model, 
the  displacement,  qUc,  z,  t;k)  of  an  isopycnal  surface  from 
equilibrium  can  be  written  as  in  Equation  {9): 
qOc,  2,  t;  k)  *=  2 


f  J 

i. 

x  j  2dn*<jt) 

-*„(*  -  V„f) 

The  function  dnja )  is  one  of  the  Jacobian  elliptic  functions 
(Press  et  al.  1992),  which  in  analogy  to  the  well-known  os¬ 
cillatory  K-dV  solution  named  “Cnoidal"  was  termed  "Dnoi* 


K©/*!3 

Figure  7.  Typical  density  profile  for  the  region  under  study. 


dal"  by  Apel  (2003).  The  function  WkJiz)r  which  gives  the 
vertical  structure  (Appendix  A),  is  an  eigenfunction  of  the 
Taylor-Goldstein  equation  with  eigenvalue  k  -  2^/^*  The 
wave  number  is  related  to  the  wavelength  at  the  trailing 
edge  of  the  packet.  A*,  as  =  2tt/Xo.  /(x,  t)  is  the  recovery 
function  that  models  the  return  to  equilibrium  over  the  semi¬ 
diurnal  tidal  cycle.  V  is  the  nonlinear  phase  velocity.  The  non¬ 
linear  parameter,  st  modulates  the  amplitude  and  the  wave¬ 
length,  thus,  defining  the  limits  of  the  packet.  The  internal 
vertical  modes  of  the  water  column  are  denoted  by  n.  The 
exponential  factor  is  an  attenuation  function  used  to  reduce 
the  amplitudes.  The  values  of  a  and  xa  are  empirically  deter¬ 
mined.  Appendix  A  provides  additional  details  regarding  the 
model. 

In  summary,  the  Dnoidal  model  reproduces  the  sharp  onset 
seen  at  the  front  of  the  packet  and,  as  the  wave  packet  pro¬ 
ceeds  in  space  and  time,  increases  packet  length  and  wave¬ 
length  while  developing  additional  oscillations.  It  also  models 
the  long-term  depression  seen  at  the  trailing  end  of  the  pack¬ 
ets.  A  study  of  the  dynamics  of  the  internal  solitons  near 
Gibraltar  using  the  Dnoidal  model  and  SAR  data  can  be 
found  in  Apel  and  Worcester  (2000).  A  comprehensive  de¬ 
scription  of  this  new  analytical  model  is  given  in  Afel  (2003)* 

Model  lnitiaIi74itioii  Using  Climatological  and  SAR 
Data 

Because  the  model  is  based  on  a  solution  to  the  K-dV  equa¬ 
tion  for  the  amplitude  and  the  Taylor-Goldstein  equation  for 
the  vertical  structure,  it  is  divided  in  two  distinct  parts.  The 
model  is  structured  so  that  the  first  module  corresponds  to 
the  Taylor-Goldstein  equation,  which  gives  the  vertical  struc¬ 
ture  (Appendix  A),  The  results  from  this  part  are  then  used 
as  input  to  a  second  module  for  the  "Dnoidal"  calculations, 
which  give  the  simulation  final  results.  In  addition  to  density 
and  current  profiles,  the  computations  of  the  first  module  re- 
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quire  an  estimation  of  the  wavelength  at  the  trailing  edge  of 
the  packet,  k*,  The  wavelength  at  the  trailing  end  of  the  pack¬ 
ets  can  be  measured  from  SAR  images  and  the  wave  number 
can  then  be  cakuJated  as  k0  ^  2^/^.  The  wave  number, 
in  turn,  is  related  to  the  amplitude,  %  to  the  nonlinear  co¬ 
efficient,  a,  and  to  the  dispersion  coefficient,  y ,  by  the  follow¬ 
ing  equation: 


This  expression  can  be  used  to  establish  consistency  between 
the  value  of  k0  and  tv  obtained  from  SAR  imagery  and  the 
parameter  values  a  and  y  given  by  the  first  section  of  the 
model. 

The  first  module  provides  the  first- mode  eigenfunction 
from  the  Taylor-Goldstein  equation,  which  gives  the  vertical 
structure  function,  the  parameter  of  no  linearity  a,  the  pa¬ 
rameter  of  dispersion  y,  and  the  linear  phase  speed,  among 
others.  TTie  second  module  in  the  model  uses  the  output  of 
the  first  section  to  carry  out  the  Dnoidal  calculations.  It  also 
requires  an  estimation  of  the  amplitude  of  the  leading  soli  ton 
pulse.  This  estimate  can  be  obtained,  using  a  two- layer  mod¬ 
el,  from  the  half-width  of  the  soli  ton  and  the  upper  and  lower 
layer  widths  in  Equation  (3).  The  half-width  of  the  soliton,  as 
explained  before,  is  another  parameter  that  can  be  extracted 
from  SAR  imagery. 

According  to  numerical  evaluations,  the  value  of  qti0  con¬ 
trols  the  number  of  oscillations  and  can  be  tuned  to  provide 
the  correct  number  of  solitons  in  the  packet.  In  general  the 
parameter  a  used  for  the  Dnoidal  calculations  has  to  be  taken 
smaller  than  the  one  calculated  by  the  models  first  module. 
This  could  be  because  during  strong  pynocline  displacements 
the  effective  depth  of  the  upper  layer  (Al)  becomes  larger, 
thus  decreasing  a  and  increasing  y,  in  Equations  (6)  and  (7). 
Afel  (2003)  suggests  that  a  more  precise  estimate  of  a  may 
be  inferred  considering  a  pynocline  depth  that,  as  the  wave 


20km 


Simulation  for  August  22nd,  1 998  N  12  4  his,  co  =  0  86357 1 


Figure  9.  Sim  uk  tion  res  u  I  ts  vs.  ub&erved  i  ate  mal  wb  ve  packet  a  fWr  one 
tidal  period-  The  area  covered  is  about  20  km  x  20  km  and  can  be  located 
in  Figure  le. 


passes,  is  displaced  by  2^,  rather  than  assuming  the  equi¬ 
librium  depth  (hi). 

In  this  simulation,  typical  density  and  current  profiles  for 
the  region  were  used  to  initialize  the  model  (Figures  7  and 
8.  The  trailing  edge  wavelength  was  estimated  at  X<,  =  320 
m,  the  half-width  of  the  leading  soliton  pulse  was  estimated 
from  SAR  imagery  at  L  -  1 13,6  m,  and  the  amplitude  was 
calculated  in  a  preceding  section  using  Equation  (3),  which 
gave  q0  —  3.3  m.  The  total  depth  was  assumed  at  60  m,  and 
the  upper  layer  depth  was  estimated  at  20  m  according  to 
previous  discussions.  The  calculations  from  the  first  module 
in  the  model  gave  a  =  “0.1072,  a  rather  high  eatimate,  and 
y  =  83,9817,  The  value  of  a  was  modified  to  a  *  -0.0397  to 


Journal  of  Coastal  Research,  Vol.  22,  No.  6,  2006 


HH 


Teixeira  tt  at 


Smujlaiton  for  August  22nd,  1998  t=24  8  ill*.  CD  =  0.663571 


Figure  10,  Simulation  results  versus  observed  internal  wave  packet  af- 
ter  two  tidal  periods.  The  area  covered  is  about  20  km  x  20  km  and  can 
be  located  in  Figure  If. 


achieve  consistency  within  the  relation  given  by  Equation 
(10).  For  these  values,  the  right  side  of  Equation  (10)  gives 
0.0008  whereas  the  left  aide  gives  0.0161,  which  are  in  fairly 
good  agreement  The  selected  value  for  a  in  combination  with 
%  =  3.3  m  provides  a  number  of  oscillations  in  harmony  with 
those  shown  in  SAR  imagery  for  the  first  two  tidal  periods 
(Figures  9  and  10).  The  linear  phase  speed  was  determined 
by  the  first  calculations  in  the  model  at  c0  =  0.8636.  Model 
results  were  monitored  for  three  tidal  periods  after  the  initial 
bore  broke  up  into  solitons. 

Validation  Studies  Using  SAR  Imagery 

The  model  was  initialized  as  described  in  the  previous  sec¬ 
tion  and  run  for  the  first  tidal  period  (12.4  h).  Results  are 
depicted  in  Figure  9,  where  they  are  compared  with  SAR  im¬ 
agery  for  the  first  soliton  packet.  The  area  covered  is  20  km 
X  20  km  and  can  also  be  seen  in  Figure  le.  There  is  consis¬ 
tency  between  observations  and  model  results  in  terms  of  the 
leading  wavelength,  number  of  soil  tons,  and  distance  trav¬ 
eled  by  the  packet.  The  number  of  soli  to  ns  seen  in  the  SAR 
image  when  looking  at  the  upper  part  of  the  picture  is  rough¬ 
ly  the  same  as  the  number  predicted  by  the  model  output. 
Note  that  there  are  usually  more  oscillations  in  a  simulated 
packet  than  the  one  appearing  in  the  SAR  image.  This  is 
because  of  SAR  resolution  problems  and  less  than  perfect 
conditions  for  image  acquisition. 

For  the  second  tidal  period  (24.8  h),  we  found  a  similar 
agreement  between  the  model  output  and  the  observed  soli¬ 
ton  train  (Figure  10).  In  this  case  the  SAR  image  only  shows 
solitons  near  the  leading  edge  and  near  the  trailing  edge  of 
the  soliton  train.  This  could  be  because  of  wind  conditions  in 
the  area  that  are  partially  masking  the  surface  manifestation 
of  the  internal  waves.  Note  that  the  surface  signature  of  the 
internal  waves  is  related  to  the  strain  rate  and,  thus,  shifted 
with  respect  to  the  internal  wave  amplitude  displacements. 
Table  3  compares  the  main  packet  characteristics  as  observed 
from  SAR  with  the  actual  model  results. 

The  slight  mismatch  in  amplitude  between  model  results 
and  observations  arises  because  the  model  was  initialised 
with  a  60  m  depth  and  with  an  amplitude  of  3.3  m,  which 
was  the  one  given  by  the  twodayer  model  equations.  After 
12.4  h,  the  Dnoidal  model  shows  that  the  amplitude  has  de¬ 
creased  to  3  nv  In  practice,  the  first  train  is  more  than  50  m 
in  depth  with  an  amplitude  of  3.7  m  whereas  the  second  train 
is  more  than  60  m  in  depth  with  an  amplitude  of  3.3  m. 

Figure  11  shows  the  evolution  of  the  sotiton  train,  as  de- 


Table  3.  Worn  packet  chamctenstics  as  observed  from  SAB  compared  with  model  results. 


Packet  Parameters 

Observations  from  SAR 

Dnoidal  Model 

First  packet  (12,4  h> 

Number  of  aohtons 

13  to  16 

20 

Distance  from  2nd  bore 

41.5  km 

42  km 

Leading  wavelength 

670  m 

700  m 

Peak  pynocline  Displacement 

3.7  m  (two-layer  model) 

3  m 

Second  Packet  (24,8  hrs) 

Number  of  solitons 

15  to  20 

30 

Distance  from  2nd  bore 

91  km 

84  km 

Leading  wavelength 

1000  m 

801m 

Peak  pynocline  Displacement 

3.3  m  (two-layer  model) 

2m 
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Figure  11,  Soli  ton  evolution  modeled  for  two  tidal  periods.  The  arrows 
correspond  to  the  centroid  of  the  packet  cj  which  moves  with  the  linear 
phase  speed  The  front  of  the  packet  moves  with  the  nonlinear  phase 
speed  c 


termined  by  the  Dnoidal  model,  at  intervals  of  6,2  h  and  dur¬ 
ing  two  tidal  periods.  The  centroid  of  the  packet  moves  with 
linear  phase  speed  c0  =  0.8636,  whereas  the  front  of  the  pack¬ 
et  moves  with  the  nonlinear  phase  speed  c  (Apel,  2003)+  This 
causes  dispersion,  something  that  is  readily  seen  in  Figure 
11. 

By  the  third  tidal  period,  the  observed  internal  wave  pack¬ 
et  shows  a  very  large  degree  of  dispersion  (Figure  lg),  which 
the  model  does  not  follow.  As  was  mentioned,  we  suspect  that 
intrusion  of  lighter  water  from  the  plume  of  the  Yangtze  Riv¬ 
er  and/or  advection  by  strong,  variable  current  fields  coupled 
with  the  highly  nonlinear  characteristics  of  this  phenomenon 
could  be  responsible  for  the  large  increase  in  wavelength  seen 
in  these  images.  An  adaptation  of  the  model  that  takes  into 
consideration  higher  nonlinearities  in  the  K-dV  equation  as 
well  as  changing  stratification  conditions  along  the  propa¬ 
gation  path  might  be  needed  to  describe  these  extreme  cases 
of  dispersion. 


CONCLUSIONS 

An  internal  wave  train  that  propagated  for  three  days  in 
the  Yellow  Sea  has  been  analyzed.  Because  part  of  the  un¬ 
derwater  topography  was  visible  through  its  sea  surface  man¬ 
ifestation  in  the  SAR  imagery,  it  provided  the  opportunity  to 
determine  with  a  very  good  degree  of  accuracy  the  origination 
point  of  the  internal  wave  train.  The  depth  profile  seen  at 
the  internal  bore  origination  region  was  used  to  initialize  a 
numerical  model  for  investigative  studies.  Several  parame¬ 
ters  of  the  solitary  internal  waves  were  derived  from  the  SAR 
data  and  used  in  coqjunetion  with  a  two-layer  model  and  the 
Dnoidal  model  to  characterize  the  dynamics  of  the  soli  tons  as 
well  as  to  infer  some  properties  of  the  water  column.  The 


model  showed  good  agreement  with  obse rv  aliens  in  the  SAR 
imagery  for  two  tidal  periods.  There  was  consistency  in  terms 
of  number  of  solitons  in  the  packet,  distance  traveled  by  the 
packet,  and  to  certain  extent,  in  the  amount  of  pynocline  dis¬ 
placement.  By  the  third  tidal  period,  the  internal  wave  pack¬ 
et  shows  a  very  large  degree  of  dispersion,  which  the  model 
does  not  follow.  An  adaptation  of  the  model  that  takes  into 
consideration  the  higher  non  linearities  in  the  K-dV  equation 
as  well  as  changing  stratification  conditions  and  ocean  cur¬ 
rents  along  the  propagation  path  might  be  needed  to  describe 
these  extreme  cases  of  dispersion.  The  large  degree  of  dis¬ 
persion  observed  within  the  internal  wave  packets  as  they 
propagated  further  southeast  was  postulated  to  be  a  possible 
consequence  of  an  increase  in  the  density  contrast  due  to  the 
seasonal  intrusion  of  the  Yangtze  River  plume  into  the  re¬ 
gion.  Following  this  hypothesis,  it  is  clear  that  sustained  ob¬ 
servations  of  internal  wave  packets  in  the  area  could  be  used 
to  remotely  monitor  regional  and  seasonal  changes  in  the  up¬ 
per  layer  relative  density  distribution.  The  study  of  further 
SAR  imagery  is  needed  to  examine  whether  the  dispersion 
extent  seen  in  these  packets  of  solitary  waves  is  associated 
with  seasonal  changes  in  stratification  conditions  or  with  ad- 
vection  by  variable  current  fields  along  the  propagation  path. 
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APPENDIX  A 

Dnoidal  Model  Equations 

The  Dnoidal  model  is  based  on  a  solution  to  the  K-dV  equa¬ 
tion  for  the  amplitude  and  the  Taylor-Goldstein  equation  for 
the  vertical  structure.  According  to  Afel  (2003),  the  displace¬ 
ment,  z>  t\ k)  of  an  isopycnal  surface  from  equilibrium 
can  be  written  in  its  final  form  as  in  the  following  equation: 
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Is  one  of  the  Jacobian  elliptic  functions,  which  de¬ 
pend  on  the  nonlinear  parameter  s ; 

Is  an  eigenvalue  of  the  Taylor-Goldstein  equation: 
k  =  2rr/X; 

Is  the  wave  number  and  is  related  to  the  wave¬ 
length  at  the  trailing  edge  of  the  packet  X*  as  - 

Zir/kfl,; 

Are  the  internal  vertical  modes  of  the  water  col¬ 
umn; 

Are  the  amplitudes  of  each  normal  mode; 

Is  the  nonlinear  phase  velocity,  which  can  be  writ¬ 
ten  as: 

v  =  e,(l  +  ■Mp-QtTlo)  (A'2) 

Is  the  nonlinear  parameter  (0  <  sa  <  1)  that  mod¬ 
ulates  the  amplitude  and  the  wavelength  thus  de¬ 
fining  the  limits  of  the  packet; 

Is  an  attenuation  fimetion  used  to  reduce  the  am¬ 
plitudes; 

It  is  empirically  determined:  a  =»  2  X  10  3  km  *; 
It  is  empirically  determined:  xa  =  1  km; 

Is  the  recovery  function  that  models  the  return  to 
equilibrium  over  the  semidiurnal  tidal  cyde: 

(A-3) 

where: 

Gives  a  scale; 

Gives  a  phase  to  the  return  to  equilibrium  profile 
at  the  trailing  end  of  the  packets;  and 
Gives  the  vertical  structure  and  is  an  eigenfunc¬ 
tion  of  the  Taylor-Goldstein  equation. 


/(x,  t)  =  - 


1  +  tanhl 


The  Vertical  Structure  Function 


The  full  Taylor-Goldstein  equation  for  linear  internal 
waves  governs  the  vertical  structure.  If  the  equation  is  gen¬ 
eralized  to  consider  a  streaming  velocity  £/0<z),  then  for  the 
eigenfunction  Wkn(z)  the  equation  becomes  that  shown  in  A-4: 


>  [  NHz)  d*U</dz* 

dz'J  +  k  [<«,  -  Jfc  t/0(2)>f  +  *K  -  A  UoW] 


0 


(A-4) 


The  TG  equation  could  be  modified  to  include  effects  of  the 
Earth's  rotation.  The  resultant  equation  is  more  complicated 
and  the  solutions  are  altered  but  mainly  at  very  low  frequen* 
cies  (Afel,  2003). 

Vertical  variations  in  the  horizontal  mean  current  velocity 
Uq(z)  and  in  the  buoyancy  frequency  N(z)  are  allowed.  For 
mode  nr  the  wave  frequency  eigenvalue  is  ov  After  imposing 
rigid-lid  boundary  conditions  at  the  surface  (z  =  9)  and  at  the 
bottom  iz  -  ~H)  this  becomes  a  second-order  Sturm-Liou- 
ville  system  (Apel,  2003), 
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